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^ ' P. Roustazadeh and M. Bottcher^ 

O 

Q ■ ABSTRACT 

The growing number of extragalactic high-energy (HE, ii^ > 100 MeV) and 
! very-high-energy (VHE, E > 100 GeV) 7-ray sources that do not belong to the 

! blazar class suggests that VHE 7-ray production may be a common property of 

^ ! most radio-loud Active Galactic Nuclei (AGN). In a previous paper, we have in- 

vestigated the signatures of Compton-supported pair cascades initiated by VHE 
7-ray absorption in monochromatic radiation fields, dominated by Lya line emis- 
Iri ! sion from the Broad Line Region. In this paper, we investigate the interaction 

of nuclear VHE 7-rays with the thermal infrared radiation field from a circum- 
nuclear dust torus. Our code follows the spatial development of the cascade in 
^ [ full 3-dimensional geometry. We provide a model fit to the broadband SED of 

^ I the dust-rich, 7-ray loud radio galaxy Cen A and show that typical blazar-like 

^ I jet parameters may be used to model the broadband SED, if one allows for an 

^ I additional cascade contribution to the Fermi 7-ray emission. 

T— I ■ 

Subject headings: galaxies: active — galaxies: jets — gamma-rays: galaxies — 
radiation mechanisms: non-thermal — relativistic processes 

1. Introduction 

Blazars are a class of radio-loud active galactic nuclei (AGNs) comprised of Flat- 
Spectrum Radio Quasars (FSRQs) and BL Lac objects. Their spectral energy distributions 
(SEDs) are characterized by non-thermal continuum spectra with a broad low-frequency 
component in the radio - UV or X-ray frequency range and a high-frequency component 
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from X-rays to 7-rays, and they often exhibit substantial variabihty across the electromag- 
netic spectrum. In the VHE 7-ray re gime, the tirae sca l e of this variability has been observed 
to be as short as just a few minutes (lAlbert et al.ll2007l : lAharonian et al.ll2007l ). While previ- 
ous generations of ground-based Atmospheric Cherenkov Telescope (ACT) facilities detected 
almost exclusively high-frequency peaked BL Lac objects (HBLs) as extragalactic sources of 
VHE 7-rays (with the notable exception of the radio galaxy M87), in recent years, a number 
of non-HBL blazars and even non-blazar radio-loud AGN have been detected by the cur- 
rent generation of ACTs. This suggests that most bl azars might be intrinsic ally emitters of 
VHE 7-rays. According to AGN unification schemes (jUrry fc Padovanilll995[ ). radio galaxies 
are the mis-aligned parent population of blazars with the less powerful PR I radio galaxies 
corresponding to BL Lac objects and PR II radio galaxies being the parent population of 
radio-loud quasars. Blazars are those objects which are viewed at a small angle with respect 
to the jet axis. If this unification scheme holds, then, by inference, also radio galaxies may 
be expected to be intrinsically emitters of VHE 7-rays within a narrow cone around the jet 



axis. 



While there is little evidence for dense radiation environments in the nuclear regions 
of BL Lac objects — in particular, HBLs — , strong line emission in Plat Spectrum Radio 
Quasars (PSRQs) as well as the occasional d etection of emission lines in the spectra of some 



BL Lac objects (e.g.. IVermeulen et al.lll995l ) indicates dense nuclear radiation fields in those 



objects. This is supported by spectral modehng of the SEDs of blazars using leptonic models 
which prefer scenarios based on external radiation fields as sources for Compt on scattering to 



produce the high-energy radiation in PSRQs, LBLs and also some IBLs ( e.g. , iGhisellini et al. 
19981 : iMadeiski et"aDll999l : iBottcher fc Bloom! I2OO0I lAcciari et aPboOsI ). If the VHE 7-ray 
emission is indeed produced in the high-radiation-density environment of the broad line 
region (BLR) and/or the d ust torus of an AGN, it is expected to be strongly attenuated 



by 77 pair production fe.g. iProtheroe fc Biermann 



1997 



Donea fc Protheroe 2003: Reimer 



20071 : iLiu et all I2OO8I : ISitarek fc Bednarekl 120081 ). lAharonian et all mm have suggested 
that such intrinsic 77 absorption may be responsible for producing the unexpectedly hard 
intrinsic (i.e., after correction for 77 absorption by the extragalactic background light) VHE 
7-ray spectra of so r ae bla zars at relatively high redshift. A similar effect has been invoked by 
Poutanen &: SternI ( 120101 ) to explain the spectral breaks in the Fermi spectra of 7-ray blazars. 
This absorption process will lead to the development of Compt o n-supported pair cascades 



in the circumnuclear environrn ent (e.g.. iBednarek fc Kirklll995l : ISitarek fc Bednarekl 12010 



Roustazadeh fc Bottcherll2010l ). 



In iRoustazadeh fc Bottchen (|2010[ ). we considered the full 3-dimensional development 
of a Compton-supported VHE 7-ray induced cascade in a monochromatic radiation field. 
This was considered as an approximation to BLR emission dominated by a single (e.g.. 
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Lya) emission line. In that work, we showed that for typical radio-loud AGN parameters 
rather small (~ /iG) magnetic fields in the central ~ 1 pc of the AGN may lead to efficient 
isotropization of the cascade emission in the Fermi energy r ange. We applied t his idea to 



fit the Fermi 7-ray emission of the radio galaxy NGC 1275 ( lAbdo et al. 



2009af ) under the 



plausible assumption that this radio galaxy would appear as a 7-ray bright blazar when 
viewed along the jet axis. 

In this paper, we pre s ent a generalization of the Monte-Carlo cascade code developed in 
Roustazadeh fc Bottcherl ( I2OIOI ) to arbitrary radiation fields. In particular, we will focus on 
thermal blackbody radiation fields, representative of the emission from a circum-nuclear dust 
torus. In Section [2] we will outline the general model setup and assumptions and describe 
the modified Monte-Carlo code that treats the full three-dimensional cascade development. 
Numerical results for generic parameters will be presented in Section |31 In Section HJ we 
will demonstrate t hat the broad-band SED of the radio galaxy Cen A, including the recent 
Fermi 7-ray data ( lAbdo et al.ll2009d ). can be modeled with plausible parameters expected 
for a mis-aligned blazar, allowing for a contribution from VHE 7-ray induced cascades in the 
Fermi energy range. We summarize in Section O 



2. Model Setup and Code Description 

Figured] illustrates the geometrical setup of our model system. We represent the primary 
VHE 7-ray emission as a mono-directional beam of 7-rays propagating along the X axis, 
described by a power-law with photon spectrum index a and a high-energy cut-off at E^^rnax- 
For the following study, we assume that the primary 7-rays interact via 77 absorption and 
pair production with an isotropic thermal blackbody radiation field within a fixed boundary, 
given by a radius -Rext, i-e., 

UeA^, r, n) = 2hiyyc^ -A H{R,^t - r) (1) 

exp(fl) - 1 

where A is a normalization factor chosen to obtain a total radiation energy density Uext (see 
Eq. [2] below), and H is the Heaviside function, H{x) = 1 if a; > and H{x) = otherwise. 
A magnetic field of order ~mG is present. Without loss of generality, we choose the y and 
z axes of our coordinate system such that the magnetic field lies in the (x,y) plane. 

The input parameters to our model simulation describing the external radiation field 
are the integral of Ue^ti^, r, fl) over all frequencies: 
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Fig. 1. — Geometry of the model setup. 



Jo 

the blackbody temperature T, and the radial extent -Rext- 



(2) 



We have used the Monte-Carlo code developed by lRoustazadeh fc Bottcherl (|2010| ). This 
code evaluates 77 absorptio n and pair production u s ing th e full analytical solution to the 
pair production spectrum of iBottcher fc Schlickeiserl (119971 ) under the assumption that the 
produced electron and positron travel along the direction of propagation of the incoming 7- 
ray. The trajectories of the particles are followed in full 3-D geometry. Compton scattering 
is evaluated using the head-on approximation, assuming that the scattered photon travels 
along the direction of motion of the electron/positron at the time of scattering. While the 
Compton energy loss to the elecron is properly accounted for, we neglect the recoil effect on 
the travel direction of the electron. In order to improve the statistics of the otherwise very 
few highest-energy photons, we introduce a statistical weight, w, inversely proportional to 
the square of the photon energy, w = Where e = 
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To save CPU time, we precalculate tables for the absorption opacity k^^, Compton 
scattering length Aic, and Compton cross section for each photon energy, electron energy 
and interaction angle before the start of the actual simulation. 

The simulation produces a photon event file, logging the energy, statistical weight, and 
travel direction of every photon that escapes the region bounded by -Rext- A separate post- 
processing routine is then used to extract angle-dependent photon spectra with arbitrary 
angular and energy binning from the log files. 
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Fig. 2. — Cascade emission at different viewing angles (/i = cos^obs)- Parameters: B = 
1 mC, Ob = 5°; Uext = 10"^ erg cm'^, i?ext = 10^^ cm, T = 1000 K, a = 2.5, E^.^ax = 5 TeV. 



Numerical Results 



For comparison with our previous stu dy on mono-energetic radiatioii fields, we conduct 
a similar parameter study as presented in ( IRoustazadeh fc Bottcherll2010l ). investigating the 
effects of parameter variations on the resulting angle-dependent photon spectra. Standard 
parameters for most simulations in our parameter study are: a magnetic field of -B = 1 mC, 
oriented at an angle of 6b = 5° with respect to the X axis {Bx = 1 mC, By = 0.1 mC); an 
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external radiation energy density of Ucxt = 10~^ erg cm~^, extended over a region of radius 
i?ext = 10^^ cm; a blackbody temperature of T = 10^ K (corresponding to a peak of the 
blackbody spectrum at a photon energy of E"^^ = 0.25 eV). The incident 7-ray spectrum has 
a photon index of a = 2.5 and extends out to -E^,max = 5 TeV. The emanating photon spectra 
for all directions have been normalized with the same normalization factor, corresponding 
to a characteristic flux level of a 7-ray bright (Fermi) blazar in the forward direction. 

Figure H] illustrates the viewing angle dependence of the cascade emission. The 77 ab- 
sorption cut-off at an energy Ec = {rrieC^Y / Eg ~ 1 TeV is very smooth in this simulation 
because of the broad thermal blackbody spectrum of the external radi ation filed. In contrast, 



the (5- function approximation for the external radiation field adopted in lRoustazadeh &: Bottcher 



( 120101 ) resulted in an artificially sharp cutoff in that work. 



At off-axis viewing angles, the observed spectrum is exclusively the result of the cas- 
cade. In the limit of low photon energies (far below the 77 absorption threshold) and 
neglecting particle escape from the cascade zone, one expects a low-frequency shape close 
to uF,, (X v^l'^ due to efficient radiati ve cooling of secondary particles injected at high en- 
ergies fiRoustazadeh fc Bottcheilboioh . However, with the typical parameters used for this 



parameter study, the assumption of negligible particle escape is not always justified. In order 
to estimate the possible suppression of the low-frequency cascade emission due to particle 
escape, we calculate the critical electron energy for which the Compton cooling time scale 
^10(7) equals the escape time scale, Tesc = -Rext/(c cos6'b). Using a characteristic thermal 
photon energy of exh = 2.8 /cT/(meC^), the resulting Compton scattered photon energy, eesc 
below which we expect to see the effects of particle escape and, hence, inefficient radiative 
cooling, is 

_ 9 X 2.8 A;r meC^ cos^^B 

corresponding to an actual energy (in GeV) of 

Ee.c ^ 2 T3 uZl cos^ Ob GeV (4) 

where T = 10^ T3 K, 

Mext = 10 ^M_5 erg cm and R = 10^*^ -Rig cm. Therefore, for our 
standard parameters, we expect the low-frequency {E < a. few GeV) to be significantly 
affected by particle escape. This explains why the low-energy photon spectra shown in 
Figure [2] are harder than z/^/^. The cascade emission is progressively suppressed at high 
energies with increasing viewing angle due to incomplete isotropization of the highest-energy 
secondary particles. This effect becomes important beyond the energy £^ic,br of Compton- 
scattered photons by secondary electrons/positrons for which the Compton-scattering length 
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\ic equals the distance travelled along the gyrational motion over an angle 9, Xicil) 
^^gyril), which is given by 



E- 



3eB 



IC,br 



4 CTt Mext 



~ 1.3 B.^uZlT^e-^ GeV. 



(5) 



where B = IQ-^ mG flRoustazadeh fc Bottcheil 120101 ) . 




E [GeV] 

Fig. 3. — The effect of a varying external radiation energy density. Parameters: B^ = 10~^ G, 
By = 1.3 X 10~^ G, 6b = 7.4°; and other parameters are the same as for Figure [2] in the 
angular bin 0.4 < /i < 0.6 

Figure |3] shows the cascade spectra for different values of the external radiation field 
energy density Uext- For large energy densities Uext ^ 10""^ erg cm~^, r^^ ^ 1 for photons 
above the pair production threshold 77 so that essentially all VHE photons will be absorbed 
and the photon fiux from the cascade becomes independent of Wgxt- 

The figure confirms our discussion concerning escape and hence inefficient radiative 
cooling of low-energy particles above (see Eq. Hj). For large values of Ucxt, the Compton 
loss time scale for all relativistic electrons producing 7-rays in the considered range, is much 
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shorter than the escape time scale. Hence, the expected uF^, cx z/^/^ shape resuhs. In the low- 
Mext case, escape affects even ultrarelativistic electrons, resulting in a substantial hardening 
of the low-energy photon spectrum. 

Figure H] illustrates the effect of a varying temperature of the external blackbody radi- 
ation field. As the temperature increases up to 1000 K the cascade flux increases because 
the 77 absorption threshold energy decreases so that an increasing fraction of 7-rays can 
be absorbed. The isotropization turnover is almost independent of T. For temperatures 
T > 1000 K the cascade flux decreases with increasing T because Uext remains fixed, leading 
to a decreasing photon number density and absorption opacity with increasing T (and, 
hence, increasing Eg). The figure also confirms our expectation (Eq. H]) that an increasing 
blackbody temperature leads to an increasing suppression of the low-frequency portion of 
the cascade emission due to particle escape. 




Fig. 4. — The effect of a varying temperature of the external blackbody radiation field. 
6b = 15° and other parameters are the same as for Figure [2]in the angular bin 0.4 < /i < 0.6 



Figures [5] and [6] illustrate the effects of varying magnetic-field parameters (strength 
and orientation). As expected, the results are essentially the same a s for c ascades in mo- 
noenergetic radiation fields investigated in iRoustazadeh fc Bottcherl ( I2OIOI ): The cascade 
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Fig. 5. — The effect of a varying magnetic field strengtli, for a fixed angle of Ob = 45° 
between jet axis and magnetic field. All other parameters are the same as for Figure |2] in 
the angular bin 0.4 < /x < 0.6 

development is extremely sensitive to the transverse magnetic field By. The limit in which 
even the highest-energy secondary particles are effectively isotropized before undergoing the 
first Compton scattering interaction, is easily reached for typical magnetic fields expected in 
the circum-nuclear environment of AGNs. 



4. Application to Cen A 



The standard AGN unification scheme (jUrry fc Padovanil 119951 ) proposes that blazars 
and radio galaxies are intrinsically identical objects viewed at different angles with respect 
to the jet axis. According to this scheme, FR I and FR 11 radio galaxies are believed to be 
the parent population of BL Lac objects and FSRQs, respectively. Hence, if most blazars, 
including LBLs and FSRQs, are intrinsically VHE 7-ray emitters potentially producing pair 
cascades in their immediate environments, the radiative signatures of these cascades might be 
observable in many radio galaxies. In fact, i^Ci^i^T provided evidence for > 100 MeV 7-ray 
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Fig. 6. — The effect of a varying magnetic field orientation, for a fixed magnetic field strength 
of i? = 1 mG. All other parameters are the same as for Figure [2] in the angular bin 0.4 < 
fi < 0.6 



emissi on from three radio galaxies (Gen A : ISreekumar et al.l (jl999l ). 3G 111: iHartman et al. 



teOOSf ). and NGG 6251: iMukheriee et al.l (120021 



high-energy 7-ray sources by Fermi (lAbdo et al.l 



The se sources have been confirmed as 
120093 ^ 



of five more radio galaxies (N GG 1275: lAbdo et al.l (l2009al ). M 87: lAbdo et al.l (l2009br ) 



20 



Obr i. alon g with the d etectio n 



3G 120, 3G 20 7, and 3G 380: lAbdo et al.l (l2010bf ll In this paper, we focus on the radio 
galaxy Gen A jAbdo et al.ll2009ch . 



The FR I Gen A is the nearest radio-loud active galaxy to Earth. It has a redshift of 
z=0. 00183 at the distance oi D = 3.7 Mpc. Recently, the Auger collaboration reported that 
the arrival directions of the highest energy cosmic rays {E > Q x 10^^ eV) observed by the 



Auger obse rvatory are correlated with nearby AGN, including Gen A (lAbraham. J., et al. 



20071 . l2008l ). This suggests that Gen A may be the dominant source of observed UHEGR 
nuclei above the GZK cut off. 



Gen A has an interesting radio structure on several size scales. The most prominent 
features are its giant radio lobes, which subtend ~ 10° on the sky, oriented primarily in 
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Fig. 7. — Fit t o the SEP of Cen A. Th e green curve is a fit to the broad-band SED using 
the model of (IBottcher fc Chiang] |2002| ). while the maroon curve is the cascade emission 
resulting from 77 absorption of the forward jet emission. The red curve is the sum of both 
contributions (viewed at an angle of 70°. 



the north-south di rection. They have been imaged a t 4.8 GHz by the Park es telescope 
(jjunkes et al.l 119931 ) and studied at up to ~ 60 GHz by iHardcastle et al.l (120091 ). The radio 
lobes are the only extrag alactic source stru cture that has so far been spatially resolved 
in GeV 7-rays by Fermi (lAbdo et al.ll2010d ). The innermost re gion of Cen A has been 
resolved with VLBI and shown to have a size of ~ 3 x 10^^ cm ( Kellermann et al.l 119971 : 



Horiuchi et al.l 120061 ). 



Meisenheimer et al. 


2007 


). 


~ 50° - 80° ( 


Tineav et al. 



infrared interferometry resolves the core size to ~ 6 x 10^^ cm 
The angle of the sub-parsec jet of Ce n A to our lin e of sight is 
1998h with a preferred value of ~ 70° flsteinlel[2009[ ). 

The K-band nuclear flux with starlight subtr acted is F(K) ~ 38 m Jy, corresponding to 
uLiy ~ 7 X 10^° erg for a distance of 3.5 Mpc (IMarconi et al. II2OOOI ). The mid-IR flux of 
~ 1.6 Jy at 11.7/im corresponds to uL^, ~ 6 x 10'^^ erg s~^. The broad band SED fr o m radi o 
to 7-rays has been fitted with a synchrotron self- Compton model by lAbdo et al.l (l2010al ). 
In their model (see Table 2 of lAbdo et al.ll2010al ). a maximum electron Lorentz factor of 
Imax = 1 X 10® was required in order to produce the observed 7-ray emission. However, 
given the assumed magnetic field of B = 6.2 G, this does not seem possible since for 7 = 10® 
electrons, the synchrotron loss time scale is shorter than their gyro-timescale, which sets 
the shortest plausible acceleration time scale. Here, we therefore present an alternative 
interpretation of the SED, based on the plausible assumption that Cen A would appear as a 
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VHE 7-ray emitting blazar in the jet direction, and cascading of VHE 7-rays on the nuclear 
infrared radiation field produces an observable off-axis 7-ray signature in the Fermi energy 
range. 



Figure [7| illustrates a broadband fit to the SED of Cen A (data from 
using the equilibrium version of the b lazar radiatio n tran sfer code of 



(120021 ). as described in more detail in lAcciari et al.l ( l2009l ). For this fit, standard blazar 



Abdo et al. 


2010a 


), 


Bottcher & Chiang 



jet parameters were adopted, but the viewing angle was chosen in accord with the observa- 
tionally inferred range. Specifically, we chose ^obs = 70°. Other model parameters include 
a bulk Lorentz factor of F = 5, a radius of the emission region of R = 1 x 10^^ cm, a 
kinetic luminosity in relativistic electrons, Lg = 9.4 x 10^^ erg s~^, a co-moving magnetic 
field of i? = 11 G, corresponding to a luminosity in the magnetic field (Poynting flux) of 
Lb = 1.1 X 10^^ erg s~^ and a magnetic-field equipartition fraction = Lb/L^ = 12, 
corresponding to a Poynting- flux dominated jet. Electrons are injected into the emission 
region at a steady rate, with a distribution characterized by low- and high-energy cutoffs at 
71 = 1.2 X 10^ and 72 = 1.0 x 10^, respectively, and a spectral index of g = 3.5. The code 
finds a self-consistent equilibrium between particle injection, radiative cooling and escape, 
from which the final photon spectrum is calculated. The resulting broadband SED fit is 
illustrated by the solid green curve in Figure [3 The flux emanating in the forward direction 
(^obs = 0°, i.e., the blazar direction) has been chosen as an input to our cascade simulation 
to evaluate the cascade emission in the nuclear infrared radiation field of Cen A observed at 
the given angle of 6'obs = 70°. 

For the cascade simulation, we assumed a blackbody temperature of 2300 K resulting 
in a peak frequency in the K-band. The external radiation field is parameterized through 
Wext = 1.5 X 10^^ erg cm~^ and i?ext = 3 x 10^^ cm. These parameters combine to an 
IR luminosity of Lr t.r = 4:7i R'^^^ cu f,^t = 5 x 10^^ erg s~^, in agreement with mid-IR flux 



observed for Cen A (IMarconi et al. II2OOOI ). The magnetic fleld is 5 = 1 mC, oriented at 
an angle of 6b = 4°. The cascade spectrum shown in Figure [7] pertains to the angular bin 
0.28 < yU < 0.38 (corresponding to 67° ^ < 73°), appropriate for the known orientation 
of Cen A and consistent with our broadband SED flt parameters. The cascade spectrum is 
shown by the maroon curve in Figure [TJ while the total observed spectrum is the solid red 
curve. The flgure illustrates that the cascade contribution in the Fermi range substantially 
improves the flt, while still allowing physically reasonable parameters for the broadband SED 
flt. 
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5. Summciry 

We investigated the signatures of Compton-supported pair cascades initiated by the 
interaction of nuclear VHE 7-rays with the thermal infrared radiation field of a circumnuclear 
dust torus in AGNs. We follow the spatial development of the cascade in full 3-dimensional 
geometry and study the dependence of the radiative output on various parameters pertaining 
to the infrared radiation field and the magnetic field in the cascade region. 

We confirm the results of our previous study of cascades in monoenergetic radiation 
fields that small (> /xG) perpendicular (to the primary VHE 7-ray beam) magnetic field 
components lead to efficient isotropization of the cascade emission out to HE 7-ray energies. 
The cascade intensity as well as the location of a high-energy turnover due to inefficient 
isotropization also depend sensitively on the energy density and temperature of the soft 
blackbody radiation field, as long as the cascade is not saturated in the sense that not all 
VHE 7-rays are absorbed. 

The shape of the low-frequency tail of the cascade emission is a result of the interplay 
between radiative cooling and escape. For environments characterized by efficient radia- 
tive cooling, the canonical uFi, oc u^^'^ spectrum results. If radiative cooling is inefficient 
compared to escape, the low-frequency cascade spectra are harder than z/^/^. 

We provide a model fit to the broadband SED of the dust-rich, 7-ray loud radio galaxy 
Cen A. We show that typical blazar-like jet parameters may be used to model the broadband 
SED, if one allows for an additional cascade contribution to the Fermi 7-ray emission due 
to 77 absorption and cascading in the thermal infrared radiation field of the prominent dust 
emission known to be present in Cen A. 

We thank J. Finke for providing the SED data points of Cen A. This work was supported 
by NASA through Fermi Guest Investigator Grants NNX09AT81G and NNX10AO49G. 
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